Fast oxidation of sulfur dioxide by hydrogen peroxide in deliquesced aerosol particles

Atmospheric sulfate aerosols have important impacts on air quality, climate, and human and ecosystem health. However, current air-quality models generally underestimate the rate of conversion of sulfur dioxide (SO2) to sulfate during severe haze pollution events, indicating that our understanding of sulfate formation chemistry is incomplete. This may arise because the air-quality models rely upon kinetics studies of SO2 oxidation conducted in dilute aqueous solutions, and not at the high solute strengths of atmospheric aerosol particles. Here, we utilize an aerosol flow reactor to perform direct investigation on the kinetics of aqueous oxidation of dissolved SO2 by hydrogen peroxide (H2O2) using pH-buffered, submicrometer, deliquesced aerosol particles at relative humidity of 73 to 90%. We find that the high solute strength of the aerosol particles significantly enhances the sulfate formation rate for the H2O2 oxidation pathway compared to the dilute solution. By taking these effects into account, our results indicate that the oxidation of SO2 by H2O2 in the liquid water present in atmospheric aerosol particles can contribute to the missing sulfate source during severe haze episodes

Accurate Measurements of Aerosol Hygroscopic Growth over a Wide Range in Relative Humidity

Using a comparative evaporation kinetics approach, we describe a new and accurate method for determining the equilibrium hygroscopic growth of aerosol droplets. The time-evolving size of an aqueous droplet, as it evaporates to a steady size and composition that is in equilibrium with the gas phase relative humidity, is used to determine the time-dependent mass flux of water, yielding information on the vapor pressure of water above the droplet surface at every instant in time. Accurate characterization of the gas phase relative humidity is provided from a control measurement of the evaporation profile of a droplet of know equilibrium properties, either a pure water droplet or a sodium chloride droplet. In combination, and by comparison with simulations that account for both the heat and mass transport governing the droplet evaporation kinetics, these measurements allow accurate retrieval of the equilibrium properties of the solution droplet (i.e., the variations with water activity in the mass fraction of solute, diameter growth factor, osmotic coefficient or number of water molecules per solute molecule). Hygroscopicity measurements can be made over a wide range in water activity (from >0.99 to, in principle, 0.9 and ∼±1% below 80% RH, and maximum uncertainties in diameter growth factor of ±0.7%. For all of the inorganic systems examined, the time-dependent data are consistent with large values of the mass accommodation (or evaporation) coefficient (>0.1)

Hygroscopic properties of aminium sulfate aerosols

Alkylaminium sulfates originate from the neutralisation reaction between short-chained amines and sulfuric acid and have been detected in atmospheric aerosol particles. Their physicochemical behaviour is less well characterised than their inorganic equivalent, ammonium sulfate, even though they play a role in atmospheric processes such as the nucleation and growth of new particles and cloud droplet formation. In this work, a comparative evaporation kinetics experimental technique using a cylindrical electrodynamic balance is applied to determine the hygroscopic properties of six short-chained alkylaminium sulfates, specifically mono-, di-, and tri-methylaminium sulfate and mono-, di-, and tri-ethyl aminium sulfate. This approach allows for the retrieval of a water-activity-dependent growth curve in less than 10 s, avoiding the uncertainties that can arise from the volatilisation of semi-volatile components. Measurements are made on particles > 5 µm in radius, avoiding the need to correct equilibrium measurements for droplet-surface curvature with assumed values of the droplet-surface tension. Variations in equilibrium solution droplet composition with varying water activity are reported over the range 0.5 to > 0.98, along with accurate parameterisations of solution density and refractive index. The uncertainties in water activities associated with the hygroscopicity measurements are typically  0.9 and  ∼  ±1 % below 0.9, with maximum uncertainties in diameter growth factors of ±0.7 %. Comparison with previously reported measurements show deviation across the entire water activity range

Hygroscopic properties of aminium sulfate aerosols

Alkylaminium sulfates originate from the neutralisation reaction between short-chained amines and sulfuric acid and have been detected in atmospheric aerosol particles. Their physicochemical behaviour is less well characterised than their inorganic equivalent, ammonium sulfate, even though they play a role in atmospheric processes such as the nucleation and growth of new particles and cloud droplet formation. In this work, a comparative evaporation kinetics experimental technique using a cylindrical electrodynamic balance is applied to determine the hygroscopic properties of six short-chained alkylaminium sulfates, specifically mono-, di-, and tri-methylaminium sulfate and mono-, di-, and tri-ethyl aminium sulfate. This approach allows for the retrieval of a water-activity-dependent growth curve in less than 10 s, avoiding the uncertainties that can arise from the volatilisation of semi-volatile components. Measurements are made on particles > 5 µm in radius, avoiding the need to correct equilibrium measurements for droplet-surface curvature with assumed values of the droplet-surface tension. Variations in equilibrium solution droplet composition with varying water activity are reported over the range 0.5 to > 0.98, along with accurate parameterisations of solution density and refractive index. The uncertainties in water activities associated with the hygroscopicity measurements are typically  0.9 and  ∼  ±1 % below 0.9, with maximum uncertainties in diameter growth factors of ±0.7 %. Comparison with previously reported measurements show deviation across the entire water activity range

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell

Hydrogen-bond structure and low-frequency dynamics of electrolyte solutions: Hydration numbers from ab Initio water reorientation dynamics and dielectric relaxation spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction

Organic component vapor pressures and hygroscopicities of aqueous aerosol measured by optical tweezers

Measurements of the hygroscopic response of aerosol and the particle-to-gas partitioning of semivolatile organic compounds are crucial for providing more accurate descriptions of the compositional and size distributions of atmospheric aerosol. Concurrent measurements of particle size and composition (inferred from refractive index) are reported here using optical tweezers to isolate and probe individual aerosol droplets over extended timeframes. The measurements are shown to allow accurate retrievals of component vapor pressures and hygroscopic response through examining correlated variations in size and composition for binary droplets containing water and a single organic component. Measurements are reported for a homologous series of dicarboxylic acids, maleic acid, citric acid, glycerol, or 1,2,6-hexanetriol. An assessment of the inherent uncertainties in such measurements when measuring only particle size is provided to confirm the value of such a correlational approach. We also show that the method of molar refraction provides an accurate characterization of the compositional dependence of the refractive index of the solutions. In this method, the density of the pure liquid solute is the largest uncertainty and must be either known or inferred from subsaturated measurements with an error of <±2.5% to discriminate between different thermodynamic treatments. (Chemical Equation Presented)

Pillared two-dimensional metal-organic frameworks based on a lower-rim acid appended calix[4]arene

Solvothermal reactions of the lower-rim functionalized diacid calix[4]arene 25,27-bis(methoxycarboxylic acid)-26,28-dihydroxy-4-tert-butylcalix[4]arene (LH₂) with Zn(NO₃)₂•6H₂O and the dipyridyl ligands 4,4/-bipyridyl (4,4/-bipy), 1,2-di(4-pyridyl)ethylene (DPE) or 4,4/-azopyridyl (4,4/-azopy) afforded a series of 2-D structures of the formulae {[Zn(4,4/-bipy)(L)]•2¼DEF}n (1), {[Zn₂(L)(DPE)]•DEF}n (2) and {[Zn(OH₂)₂(L)(4,4/-azopy)]•DEF}n (3) (DEF = diethylformamide)

Modelling the hygroscopic growth factors of aerosol material containing a large water-soluble organic fraction, collected at the Storm Peak Laboratory

The compositions of six aggregated aerosol samples from the Storm Peak site have been comprehensively analysed (Hallar et al., 2013), focusing particularly on the large water-extractable organic fraction which consists of both high molecular weight organic compounds and a range of acids and sugar-alcohols. The contribution of the soluble organic fraction of atmospheric aerosols to their hygroscopicity is hard to quantify, largely because of the lack of a detailed knowledge of both composition and the thermodynamic properties of the functionally complex compounds and structures the fraction contains. In this work we: (i) develop a means of predicting the relative solubility of the compounds in the water-extractable organic material from the Storm Peak site, based upon what is known about their chemical composition; (ii) derive the probable soluble organic fraction from comparisons of model predictions with the measured hygroscopicity; (iii) test a model of the water uptake of the total aerosol (inorganic plus total water-extractable organic compounds). Using a novel UNIFAC-based method, different assignments of functional groups to the high molecular weight water soluble organic compounds (WSOC) were explored, together with their effects on calculated hygroscopic growth factors, constrained by the known molecular formulae and the double bond equivalents associated with each molecule. The possible group compositions were compared with the results of ultrahigh resolution mass spectrometry measurements of the organic material, which suggest large numbers of alcohol (–OH) and acid (–COOH) groups. A hygroscopicity index (HI) was developed. The measured hygroscopic growth is found to be consistent with a dissolution of the WSOC material that varies approximately linearly with RH, such that the dissolved fraction is about 0.45–0.85 at 90% relative humidity when ordering by HI, depending on the assumptions made. This relationship, if it also applies to other types of organic aerosol material, provides a simple approach to calculating both water uptake and CCN activity (and the κ parameter for hygroscopic growth). The hygroscopicity of the total aerosol was modelled using a modified Zdanovskii-Stokes-Robinson approach as the sum of that of the three analysed fractions: inorganic ions (predicted), individual organic acids and “sugar alcohols” (predicted), and the high molecular weight WSOC fraction (measured). The calculated growth factors broadly agree with the measurements, and validate the approach taken. The insights into the dissolution of the organic material seem likely to apply to other largely biogenic aerosols from similar remote locations

Chemical speciation models based upon the Pitzer activity coefficient equations, including the propagation of uncertainties. II. Tris buffers in artificial seawater at 25 °C, and an assessment of the seawater ‘Total’ pH scale

The substance Tris (or THAM, 2-amino-2-hydroxymethyl-1,3-propanediol, CAS 77–86-1), and its protonated form TrisH+, is used in the preparation of pH buffer solutions for applications in seawater chemistry. The development of an acid-base chemical speciation model of buffer solutions containing Tris, TrisH+, and the major ions of seawater is desirable so that: (i) the effects of changes in the composition of the medium on pH can be calculated; (ii) pH on the free (a measure of [H+]) and total (a measure of ([H+] + [HSO4−])) scales can be interconverted; (iii) approximations inherent in the definition of the total pH scale can be quantified; (iv) electrode pairs such as H+/Cl− and H+/Na+ can more easily be calibrated for the measurement of pH. As a first step towards these goals we have extended the Pitzer-based speciation model of Waters and Millero (Mar. Chem. 149, 8–22, 2013) for artificial seawater to include Tris and TrisH+, at 25 °C. Estimates of the variances and covariances of the additional interaction parameters were obtained by Monte Carlo simulation. This enables the total uncertainty of any model-calculated quantity (e.g., pH, speciation) to be estimated, as well as the individual contributions of all interaction parameters and equilibrium constants. This is important for model development, because it allows the key interactions to be identified. The model was tested against measured EMFs of cells containing Tris buffer in artificial seawater at 25 °C, and the mean deviation was found to be 0.13 ± 0.070 mV for salinities 20 to 40. Total variances for calculated electromotive forces of the buffer solutions are dominated by contributions from just a few interaction parameters, making it likely that the model can readily be improved. The model was used to quantify the difference between various definitions of total pH and –log10([H+] + [HSO4−]) in Tris buffer solutions at 25 °C, for the first time (item (iii) above). The results suggest that the total pH scale can readily be extended to low salinities using the established approach for substituting TrisH+ for Na+ in the buffer solutions, especially if the speciation model is used to quantify the effect on pH of the substitution. The relationships between electromotive force (EMF), and pH on the total scale, with buffer molality in artificial seawater at constant salinity are shown to be linear above about 0.01 to 0.02 mol kg−1 buffer molality. The pH of Tris buffers containing ratios of TrisH+ to Tris that vary from unity can be calculated very simply. Technical aspects of the total pH scale, such as the extrapolation of pH to zero buffer (at constant salinity), are clarified. Recommendations are made for further work to extend the model to the temperature range 0–45 °C, and improve accuracy, so that requirements (i) to (iv) above can be fully met